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LONG-TERM  BIOEFFECTS  OF  435-MHZ  RADIOFREQUENCY  RADIATION 
ON  SELECTED  BLOOD-BORNE  ENDPOINTS  IN  CANNULATED  RATS 
VOLUME  2.  Plasma  ACTH  and  Plasma  Corticosterone 

I.  INTRODUCTION 

During  the  past  50  years*  the  United  States  has  witnessed  a  period  of 
explosive  growth  In  the  radar  and  communications  fields.  This  growth  has 
increased  the  demand  for  available  bandwidth*  and  this  has  pushed  radar  and 
communications  frequencies  Into  higher  and  higher  ranges.  Higher  frequency 
ranges  have  permitted  faster  data  transmission  rates  and  reduced  Intersystem 
electromagnetic  Interference.  However*  these  advances  have  come  at  the  expense 
of  altering  the  planet’s  radiofrequency  radiation  (RFR)  environment.  Until  the 
advent  of  advanced  radar  and  communications*  cosmic  rays  and  background 
radiation  were  the  primary  sources  of  the  Earth’s  electromagnetic  environment. 
Radar  and  communications  transmissions  have  since  Increased  the  electromagnetic 
background  or  ambient  radiation  at  the  planet’s  surface  by  several  orders  of 
magnitude.  At  this  time*  the  biological  effects  of  exposure  to  this  omnipresent 
electromagnetic  environment  are  not  well  understood*  despite  studies  conducted 
over  the  past  several  decades. 

This  report  presents  results  of  plasma  adrenocorticotropic  hormone  (ACTH) 
and  plasma  corticosterone  assays  of  blood  samples  drawn  from  a  large  population 
of  male  Sprague-Dawley  rats  exposed  to  a  1.0  mW/cm^»  435-MHz  pulsed-wave  (1.0-ys 
pulse  width,  1-kHz  pulse  rate)  RFR  environment  for  a  6-month  duration.  The 
exposure  group  consisted  of  100  cannulated  rats  housed  In  Plexiglas  cages 
arrayed  on  the  tiers  of  a  stacked*  parallel-plate  circular  waveguide. 
Engineering  aspects  of  this  waveguide  and  the  exposure  environment  It  generated 
have  been  previously  reported  [13.  The  sham-exposure  group  consisted  of  100 
cannulated  rats  housed  in  an  Identical,  but  unenergized*  collocated  facility. 
Results  reporting  plasma  prolactin  concentrations  In  these  same  animals  are  in 
Volume  3  of  this  series  [23. 

The  pituitary  gland  releases  a  number  of  hormones  that  regulate 
physiological  functions  in  the  body.  One  of  these  hormones* 
adrenocortlcotrophlc  hormone  (ACTH)*  stimulates  the  adrenal  cortex  to  release 
corticosterone  and  other  cortical  hormones,  such  as  cortisol.  In  rats  and  mice, 
only  corticosterone  Is  secreted.  The  primary  physiologic  role  of  ACTH  is  to 


1 


stimulate  the  secretion  and  synthesis  of  corticosterone  (and  of  cortisol  In  man 
and  In  many  animal  species  except  the  rat)  by  the  adrenal  cortex.  ACTH  also 
acts  as  a  trophic  substance  maintaining  the  size  and  blood  flow  to  the  adrenal 
cortex  and  exhibits  adrenal  physiologic  effects  on  cyclic-adenosine 
monophosphate  (c-AMP)  mediated  systems.  Both  ACTH  and  corticosterone  are 
readily  affected  by  stress.  In  all  species  observed  so  far  and  In  all 
experimental  situations*  a  stressful  condition  leads  to  an  increased  secretion 
of  ACTH  [3]  and  of  corticosterone  [4,53.  The  increase  In  stress  hormones 
corresponds  to  the  stimulus*  being  a  graded  function  of  stimulus  Intensity  and 
duration  [6],  Depending  on  intensity  (and  duration)  of  stress*  plasma 
concentration  of  the  hormones  Is  Increased  several-fold.  This  fact  permits  the 
Increase  of  plasma  ACTH  and  corticosterone  concentration  to  be  used  to  measure 
stress  Intensity.  In  rats*  levels  of  plasma  corticosterone  vary  from  6-8  ug/100 
mL  in  unstressed  animals  to  80  pg/100  mL  in  stressed  animals*  a  tenfold 
Increase.  This  Increase  Indicates  how  sensitive  the  hypothalamic-pituitary- 
adrenal  axis  Is  to  physiological  and  psychological  Insults  [7], 

Anesthesia,  exercise  [8],  Immobilization  [9],  withdrawal  of  large  volumes 

of  blood  [10],  exposure  to  a  new,  unfamiliar  cage  or  room  [11,12]*  noise, 

hypoxia  [13*14],  a  decreased  Pan  *  handling  of  the  animals  [15],  cold  exposure 
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[16],  and  many  other  environmental  factors  lead  to  Increases  In  plasma 
concentration  of  ACTH  and  corticosterone  and,  therefore*  had  to  be  avoided  In 
this  study.  Even  smoking  In  the  animal  housing  room  Increases  rat  plasma  ACTH 
and  plasma  corticosterone  levels  [17],  Clearly*  both  neurogenic  (emotional)  and 
systemic  (somatic)  stimuli  are  effective  In  evoking  Increased  ACTH  and 
corticosterone  secretion  In  animals  and  In  man. 

Circadian  rhythm  Is  a  known  factor  affecting  plasma  hormone  concentrations 
In  rats.  Both  plasma  ACTH  and  plasma  corticosterone  follow  circadian  rhythm, 
rising  each  day  during  the  evening  hours  and  decreasing  to  the  lowest  level 
between  9  AM  and  1  PM  [18,19]. 

Some  researchers  have  suggested  that  individual  housing  might  lead  to 
Increases  In  plasma  corticosterone  In  rats  (called  "Isolation-stress  syndrome," 
[20*21,22]).  However*  most  Investigators  found  the  difference  between  slngle- 
or  group-housed  rats  to  be  very  small  or  nonexistent.  Results  of  experiments 
during  this  study  also  Indicated  that  values  of  resting  plasma  ACTH  and  plasma 
corticosterone  were  similar  In  both  single-  and  In  group-housed  rats. 
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The  ACTH  and  corticosterone  Increase  In  response  to  stress  can  be 
quantified.  The  degree  of  plasma  ACTH  or  plasma  corticosterone  Increase  Is 
related  to  the  type  of  stress  to  which  the  animal  Is  exposed  [11]. 
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II.  MATERIALS  AND  ►CTHODS 


For  this  study*  the  concentrations  of  both  plasma  ACTH  and  plasma 
corticosterone  were  chosen  as  sensitive  indicators  of  possible  environmental 
stresses  Induced  by  RFR.  To  detect  and  quantitatively  evaluate  possible 
Increases  In  plasma  ACTH  and  plasma  corticosterone  levels  Induced  by  RFR,  blood 
was  sampled  and  assayed  from  86  exposed  and  65  sham-exposed  animals  (in  the  case 
of  ACTH);  87  exposed  and  65  sham-exposed  animals  (in  the  case  of 
corticosterone).  The  animals  assayed  for  ACTH  concentration  correspond  to  the 
animals  assayed  for  corticosterone  concentration.  Analysis  of  the  data  obtained 
from  the  blood  sample  assays  determined  whether  there  were  any  RFR-induced 
changes  In  plasma  ACTH  and  plasma  corticosterone  concentrations. 

Animals.  The  rat  represents  a  comparatively  Inexpensive  and  homogeneous 
population.  For  this  reason.  It  Is  often  desirable  to  use  this  species  as  the 
animal  model  In  physiologic  studies. 

We  used  male  Sprague-Dawley  rats  In  this  study.  All  experimental  animals 
were  obtained  from  the  same  building  and  room  at  CAMM  Research  Labs,  Wayne,  New 
Jersey.  The  animals,  weighing  approximately  60  g,  were  delivered  to  Emory 
University  where  they  were  caged  singly  and  given  water  and  food  (Purina  Rat 
Chow)  ad  1  Ibltum.  Temperature  In  the  animal  rooms  was  maintained  at  24  +  I  °C 
and  the  photoperiod  was  12  hours/12  hours,  with  the  lighted  phase  occurring 
between  8  AM  and  8  PM. 

Experimental  Facility.  The  Georgia  Tech  Research  Institute’s 
Radiofrequency  Radiation  Facility  [23]  consisted  of  8  collocated  rooms  on  the 
basement  floor  of  the  Baker  Building  on  the  main  campus.  These  eight  rooms 
provided  a  closed,  complete  facility  for  long-term  bioeffects  studies  Involving 
rodents. 

Two  Identical,  collocated  rooms  In  the  Facility  housed  the  100  exposure  and 
100  sham-exposure  animals.  Each  room  contained  a  stack  of  circular,  parallel- 
plate  waveguides  fed  by  a  slotted-cyl Inder  antenna  system  for  radiating  the 
animals.  The  stacks  of  parallel  waveguides  consisted  of  five  3.6-m  (12  ft)  dia. 
plates  that  made  up  4  sets  of  circular  waveguides.  Twenty-five  Individually 
housed  rats  were  positioned  around  the  circumference  of  each  waveguide  set.  The 
walls  of  both  rooms  were  lined  with  anecholc  absorbing  material  and  shielded 
with  aluminum  foil  to  prevent  excessive  microwave  leakage  radiation. 
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The  circular*  parallel-plate  waveguide  assembly  provided  a  1.0  mW/cnr 
exposure  field  around  the  circumference  of  the  plates.  The  45.7-cm  (18  In.) 
plate  separation  distance  permitted  propagation  of  a  TE^q  mode  wave  with 
horizontal  polarization.  The  power  density  displayed  a  cosine-squared 
dependency  between  the  plates*  with  the  maximum  power  density  occurring  midway 
between  each  set  of  plates.  This  arrangement  positioned  the  electric  field 
vector  parallel  to  the  rat's  longitudinal  axis*  thereby  maximizing  the  coupling 
between  the  electric  field  and  the  rat. 

A  si otted-cyl  1  nder  antenna  with  the  proper  diameter*  thickness*  slot 
length*  and  slot  width  dimensions  fed  the  stack  of  circular  waveguides  In  a 
manner  that  provided  an  essentially  constant  electric  field  Intensity  In  the 
azimuth  plane. 

Cages.  The  cages  were  constructed  of  Plexiglas  to  facilitate  visual 
observation  of  the  rats.  Each  cage  was  22.9-cm  (9  in.)  long  by  12.7-cm  (5  in.) 
wide  by  17.8-cm  (7  in.)  tall.  These  dimensions  complied  with  dimensions 
recommended  by  the  National  Institutes  of  Health  for  long-term  housing  of  rats 
[24].  The  food  hopper  and  water  bottle  were  placed  on  the  distal  side  of  the 
cage  to  minimize  their  Interaction  with  the  exposure  field.  The  glass  floor 
rods  In  the  cage  were  oriented  perpendicular  to  the  cage's  long  axis  to  Induce 
the  rats  to  preferentially  align  themselves  parallel  to  the  electric  field 
vector.  The  sipper  tubes  of  the  water  bottles  were  made  of  glass  to  be 
nonperturbing  In  the  field.  Evaluations  of  the  cages  conducted  In  the  circular* 
parallel-plate  waveguide  assembly  showed  field  scattering  from  the  Plexiglas  and 
water  to  be  below  the  range  of  detection. 

The  RFR  Facility  contained  a  data  acquisition  system  for  storing  and 
processing  experimental  data*  an  electronic  balance  for  weighing  the  rats  during 
the  study*  and  rooms  for  transmitter  operation*  blood  sampling*  cage  washing* 
and  materials  storage. 

Noise  Increases  the  concentration  of  both  plasma  ACTH  and  plasma 
corticosterone  In  rats  (Fig.  1).  To  avoid  the  possible  effects  of  noise  during 
this  study*  the  entire  RFR  Facility  was  kept  locked  to  avoid  unauthorized  entry. 
Only  the  animal  caretaker  and  the  technician  who  sampled  blood  from  the  animals 
were  permitted  uncontrolled  entry  to  the  Facility. 

Cannulatlon.  To  detect  and  quantitatively  evaluate  Increases  in  plasma 
ACTH  and  plasma  corticosterone*  the  resting  levels  of  these  hormones  first  had 
to  be  determined.  It  became  Immediately  clear  that*  to  obtain  the  real  resting 
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Figure  1.  Effect  of  loud  on-off  noise  on  plasma  corticosterone  In 
cannulated  rats  (N  =  10). 


values  of  these  two  hormones  In  undisturbed  animals#  many  routine  techniques  for 
handling  the  animals  and  for  sampling  the  blood  could  not  be  used.  For  example, 
guillotine  blood  sampling  techniques  commonly  employed  In  many  endocrinological 
studies  were  not  suitable  for  this  study.  To  use  each  animal  as  Its  own 
control,  arterial  blood  was  sampled  by  means  of  chronically  implanted  aortic 
cannulas  [25,26].  This  simple.  Inexpensive  technique  permitted  remote  stress- 
free  blood  sampling  In  conscious,  unrestrained,  and  resting  rats.  Arterial  blood 
drawn  from  the  chronically  Implanted  cannulas  was  assayed  for  plasma  ACTH  and 
plasma  corticosterone. 

The  Idea  of  sampling  venous  blood  from  the  animals  was  abandoned.  In 
venous  blood  vessels,  the  flow  regime  is  laminar  with  blood  flowing  In  discrete 
layers.  The  layers  of  blood  In  the  middle  of  the  vessels  travel  much  faster 
than  those  close  to  the  vessel  walls.  The  most  Important  consideration, 
however,  was  that  blood  layers  do  not  mix  In  venous  blood  vessels.  Thus,  a 
sample  of  venous  blood,  withdrawn  with  a  needle  or  a  cannula,  might  represent 
the  blood  returning  from  one  part  of  the  body  or  the  other,  from  a  single  organ 
or  muscle,  or  from  any  one  of  the  endocrine  glands.  For  this  reason,  we  decided 
to  sample  arterial  blood,  which  Is  always  fully  mixed.  The  mixing  occurs  In  the 
left  ventricle  of  the  heart  and  In  early  parts  of  the  aorta.  Only 
physiologically  minute  amounts  of  arterial  blood  (up  to  0.3  mL)  were  withdrawn 
from  resting  rats  approximately  once  every  2  weeks. 
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PE-10  arterial  cannulas  were  used  In  this  study.  Larger  PE-50  cannulas 
were  unsuitable  because  they  could  develop  large  blood  clots  If  not  drained 
frequently.  Large  cannulas  require  multiple  flushing  to  remain  patent#  but 
flushing  might  Induce  multiple  strokes  In  the  animals.  Use  of  1000  unlts/mL 
heparin  In  the  cannula  tip  allowed  spontaneous  patency  of  PE-50  cannulas  to 
persist  only  up  to  2  days.  Chronic  cannulation  of  the  aorta  with  a  PE-10 
cannula  was  preferable  to  cannulation  of  other  arterial  blood  vessels.  The  need 
for  properly  mixed  blood  precluded  cannulation  of  either  the  jugular  vein  or  the 
Inferior  vena  cava.  Cannulation  of  the  abdominal  aorta  provided  long-term 
functional  cannulas  [27],  Ljt  the  cannulation  procedure  was  lengthy  (20-30  min) 
and  required  temporary  dislocation  of  the  Intestinal  system.  The  abdominal 
aortic  cannula  had  a  much  larger  dead  space  than  the  aortic  cannula. 
Cannulation  of  the  aorta  through  the  left  carotid  artery,  on  the  other  hand, 
required  an  incision  of  1-1.5  cm  that  neither  penetrated  body  walls  nor  entered 
the  abdominal  cavity.  Further,  this  cannulation  could  be  completed  In  about  8 
min. 

The  carotid  artery  of  the  animal  was  cannulated  8  to  10  days  before  the 
animals  entered  the  study.  The  surgery  was  done  using  ketamlne-xylazlne 
anesthesia  (1:1  mixture;  ketamine  100  mg/mL,  xylazlne  20  mg/mL,  Intramuscular 
0.1  mL/100  g  of  body  weight).  The  catheter  was  filled  with  slightly  heparinized 
saline*  and  the  distal  end  was  sealed  with  a  nylon  plug.  Stress  hormone  levels 
returned  to  the  basal  values  approximately  3  days  after  Implantation  of  the 
chronic  arterial  cannulas.  The  first  blood  sampling  occurrad  10  days  after 
aortic  cannulation. 

Blood  Sampling.  Restraint  and  handling  Increase  stress  hormone  levels  In 
rats  (Fig.  2).  This  Increase  will  persist  20  to  30  min  after  the  Initial 
stimulus,  and  Is  apparent  even  if  the  stimulus  Is  removed  Immediately  [7], 
However,  the  animals  had  to  be  handled  upon  removal  from  their  exposure  cage  and 
placement  In  the  "sampling  box"  In  preparation  for  blood  withdrawal.  To  avoid 
the  undesired  effects  of  handling  and  stress  on  hormone  levels,  blood  from  the 
aortic  cannula  was  sampled  30  min  after  the  animal  was  placed  In  the  sampling 
box.  This  procedure  permitted  the  altered  plasma  ACTH  and  plasma  corticosterone 
levels  sufficient  time  to  return  to  their  basal  (resting)  values.  Each  animal 


*0.5  cm3 heparin  sodium  (from  beef  lung),  1000  unlts/mL  per  30  cm3 saline. 
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Figure  2.  Light  restraint  (placement  of  rats  In  narrow  Plexiglas 
cages)  Increases  plasma  corticosterone  level  of  cannulated 
rats  (N  =  10). 


was  preconditioned  for  the  sampling  box  through  a  regime  of  several  30-m1n*-long 
experiments  conducted  during  a  1-week  period  before  entering  the  study. 

After  acclimating  for  30  min  In  the  sampling  box,  the  rat's  cannula  was 
positioned  through  the  slot  In  the  top  of  the  box  (Fig.  3).  The  heparinized 
saline  was  then  removed  from  the  cannula,  and  a  0.3-mL  blood  sample  was  taken 
from  the  resting  rat.  The  withdrawal  of  larger  amounts  of  blood  from  the 
cannulated  rats  led  to  Increased  plasma  corticosterone  levels  (Fig.  4).  Using  a 
sterile  l-cm3tubercul In  syringe  fitted  with  a  30-ga  needle,  the  blood  sample  was 
taken  from  the  cannula.  The  syringe  and  the  needle  were  rinsed  with 
ethylenedlamlnetetracetate  (EDTA)  before  sampling.  The  blood  sample  was  placed 
In  an  EDTA-treated  0.3-mL  capillary  blood  collection  container  (Walter  Sarstedt 
Co.,  Princeton,  New  Jersey),  shaken,  and  then  placed  on  Ice.  The  blood  sampling 
procedure  required  about  2  min  for  each  rat. 

Anesthesia  was  not  used  during  blood  sampling  because  anesthetic  agents 
affect  plasma  stress  hormone  concentrations.  In  this  study,  an  increase  In 
plasma  corticosterone  was  observed  after  barbiturate  anesthesia  (Fig.  5). 

Plasma  ACTH  and  corticosterone  both  follow  circadian  rhythm  In  the  rat 
(Fig.  6).  To  avoid  any  possible  effect  of  the  circadian  rhythm,  blood  sampling 
occurred  between  9  AM  and  1  PM  only,  the  period  when  both  stress  hormones  were 
at  their  lowest  (true  resting)  level. 
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Figure  4.  At  -20  min  10  cannulated  rats  were  removed  from  their  home 
cages  and  placed  In  sampling  boxes.  At  this  moment  a  first 
0.2  mL  sample  of  blood  was  withdrawn.  At  0  time  (20  min 
later)  another  0.2  mL  sample  of  blood  was  withdrawn  for 
plasma  corticosterone  determination  and  the  animals  were 
bled  1  mL  (blood  withdrawn  from  cannula). 
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PLASMA  CORTICOSTERONE  C^g/IOOml)  PLASMA  CORTICOSTERONE  j^g/IOOml) 


ure  5.  Plasma  corticosterone  In  resting  cannulated  rats  and  In  the 
same  animals  10  min  after  administration  of  light  Nembutal 
anesthesia  (2  groups*  10  animals  each). 


Figure  6.  Circadian  rhythm  of  plasma  corticosterone  In  cannulated  rats 
(N=10) . 
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Blood  Sampling  Schedule.  F1gur3  7  shows  the  sampling  schedule  designed  for 
the  experiment.  Note  that  the  200  rats  were  Introduced  into  the  study  In  4 
groups  of  50  animals  each.  The  groups  entered  in  a  staggered  manner  to 
facilitate  the  process  of  logging  in  and  establishing  the  new  animals.  Each 
group  contained  25  exposure  and  25  sham-exposure  animals.  Of  the  25  exposure 
(or  sham-exposure)  animals*  20  were  sampled  for  plasma  stress  hormones*  while 
the  remaining  5  were  used  for  hematology  studies. 

The  sampling  duration  was  36  weeks  long*  Including  a  6-week  preexposure 
adaptation  period*  a  24-week  exposure  period*  and  a  6-week  postexposure  period. 
Allowing  for  group  staggering,  the  experiment  duration  was  42  weeks  long  (since 
Introduction  of  the  4  groups  was  staggered  In  2-week  Intervals).  Plasma  ACTH 
and  plasma  corticosterone  were  sampled  for  all  time  periods  marked  (A)  In  Figure 
7.  Therefore,  each  animal  should  have  been  sampled  fcr  both  plasma  ACTH  and 
plasma  corticosterone  at  weeks  -6,  -3*  0,  3,  6,  ...»  27.  This  schedule  was 
rather  rigorous,  and  therefore  could  tolerate  slight  fluctuations  In  protocol 
without  111  effects. 
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ACTH  and  Corticosterone  Determinations.  Radioimmunoassays  are  rapid, 
sensitive,  specific,  reliable,  and  require  a  minimum  quantity  of  blood.  The 
last  feature  was  very  Important,  especially  since  repetitive  sampling  was 
required  and  small  laboratory  animals  were  used. 

The  first  radioimmunoassay  was  developed  and  Introduced  In  1968  by  Berson 
and  Yalow  [283.  The  technique  was  based  on  the  principle  of  Isotopic  dilution 
In  the  presence  of  specific  antibodies.  Although  many  variants  exist,  classical 
radioimmunoassays  of  hormones  depend  upon  the  competition  between  radiolabeled 
and  unlabeled  analyte  (antigen  or  antigenic  substance)  for  a  limited  number  of 
binding  sites.  Provided  the  binding  affinities  for  the  labeled  and  unlabeled 
antigens  are  Identical#  the  quantities  of  each  antigen  bound  are  directly 
proportional  to  their  relative  concentrations.  The  total  concentration  of 
labeled  antigen  Is  the  same  In  each  Incubation  mixture  of  reagents.  These 
reagents  are  readily  available  from  commercial  sources.  Reagents  for  this  study 
were  obtained  from  Immunonuclear  Corporation.  Stillwater.  Maine.  With  this 
assay,  normal  ACTH  concentration  could  be  measured  In  unextracted  plasma.  The 
lower  limit  of  sensitivity  was  determined  by  the  extent  to  which  Incubation 
damage  during  the  assay  Incubation  period  could  be  controlled,  since  this 
affected  the  amount  of  plasma  that  the  Incubation  mixture  could  tolerate. 

Plasma  ACTH  Assay.  Following  Instructions  provided  by  Immunonuclear 
Corporation,  the  plasma  ACTH  assay  procedure  was  as  follows: 

1.  Set  up  and  label  12  x  75  mm  disposable  glass  tubes  In  duplicate  for  each 

standard.  Including  the  zero,  each  control,  and  each  sample. 

2.  Place  the  pack  of  tubes  on  crushed  Ice. 

3.  Add  100  pl_  of  each  standard  Including  zero,  control  or  experimental  sample 
to  the  appropriate  tubes. 

4.  Add  200  pL  of  rabbit  antl-ACTH  serum  In  all  tubes  except  total  count  tubes. 

5.  Add  200  pL  of  125I-ACTH  to  all  tubes. 

6.  Mix  and  Incubate  for  24  h  at  2-8  °C. 

7.  Add  500  pL  goat  anti-rabbit  precipitating  complex  to  each  tube. 

8.  Mix  and  Incubate  for  15-25  min  at  20-25  °C. 

9.  Centrifuge  for  20  min  at  760  x  g  at  20-25  °C. 

10.  Immediately  decant  the  supernate  from  all  tubes  except  the  total  count 
tubes  by  Inverting  them  for  a  maximum  of  2  min. 
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11.  In  a  gamma  scintillation  counter#  count  the  precipitate  of  each  tube  and 
the  total  count  tubes  for  a  sufficient  time  to  achieve  statistical 
accuracy. 

12.  Run  samples  with  0#  28#  63#  120#  240#  and  600  pg/mL  standards  and  three 
commercial  controls  at  the  beginning  and  at  the  end  of  the  assay.  In 
addition#  Include  three  pooled  rat  plasma  samples.  Standard  curves  are 
plotted  and  values  determined  by  a  computer  program. 


Plasma  Corticosterone  Assay.  The  assay  for  corticosterone  was  a  direct 
radioimmunoassay  for  corticosterone  which  required  25  yL  of  rat  plasma.  The 
method  was  based  on  a  procedure  for  cortisol  determination  In  man  described  by 
Donohue  and  Sgoutas  [29].  After  heat  Inactivation  of  corticosterone  binding 
proteins  In  rat  plasma#  plasma  corticosterone  was  assayed  with  specific  anti¬ 
sera  raised  against  cortlcosterone-21-heml succinate.  While  this  antibody  had  a 
significant  cross  reactivity  with  cortisol#  prior  separation  of  corticosterone 
was  not  necessary  since  the  rat  does  not  secrete  cortisol. 

The  procedural  steps  Included; 

A.  Preparation  of  reagents 

1.  An  assay  buffer  which  was  0.05  M  Trls#  0.1  M  NaCl#  and  0.1%  BSA#  pH 

8.0. 

2.  A  boiling  buffer  which  was  0.05  M  Trls  and  0.1  M  NaCl#  pH  8.O.- 

3.  (^H)  corticosterone  (0.25  m  C1»  spec.  act.  80-100  Cl/mmol;  New  England 
Nuclear  Corp.)  which  was  diluted  to  25  mL  with  absolute  ethanol.  For 
the  working  solution#  0.1  mL  was  diluted  to  10  mL  with  assay  buffer. 

4.  An  antiserum  which  was  purchased  from  Miles  Laboratories.  Each  vial 
was  reconstituted  with  5  mL  assay  buffer.  The  volume  of  the 
reconstituted  antiserum  was  diluted  with  an  equal  volume  of  assay 
buffer  for  the  working  solution. 

5.  Dextran-coated  charcoal. 

6.  Preparation  of  corticosterone  standard  by  dissolving  10  mg  of 
corticosterone  in  500  mL  ethanol. 


B.  Procedure 

1.  Standards.  To  a  series  of  tubes  labeled  A  through  G#  In  duplicate# 
add  0.1  mL  of  corticosterone  standard  to  tube  A#  evaporate  ethanol 
under  a  stream  of  nitrogen  or  air#  then  add  1.25  mL  of  assay  buffer. 
Vortex  well.  To  the  tubes  labeled  B  through  G#  add  0.5  mL  of  assay 
buffer.  Add  0.5  mL  of  the  contents  of  A  to  tube  B.  Vortex  and 
transfer  0.5  mL  from  B  to  the  next  tube  and  so  on  until  the  following 
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serial  dilutions  of  corticosterone  are  prepared:  160*  80 #  40,  20#  10# 
5  and  2.5  pg/100  mL. 


2.  Prepare  specimens  and  controls  for  assays  as  follows:  Add  0.025  mL  of 

sample  and  1.0  mL  of  boiling  buffer  In  duplicate  to  appropriately 

labeled  tubes. 

3.  Cap  the  standard  and  sample  tubes  loosely  and  place  them  In  a 

vigorously  boiling  water  bath  for  30  min. 

4.  Label  a  complete  set  of  tubes  Including  tubes  for  standards#  samples# 

total  counts#  and  zero  point  or  maximum  binding#  and  counts  not 
absorbed  by  the  charcoal  In  the  absence  of  antibody  or  background. 

5.  Add  0.05  mL  of  each  standard#  sample#  or  control  to  the  appropriate 

tubes.  Add  0.1  mL  of  radioactive  corticosterone  working  solution  to 
all  tubes.  Add  0.1  mL  of  the  antiserum  working  solution  Into  all 
assay  tubes  except  background  and  total  count  tubes. 

6.  Vortex  tubes  and  place  them  In  a  37  °C  water  bath  for  60  min. 

7.  Place  tubes  In  an  Ice  bath  for  15  min. 

8.  Allow  Dextran-coated  charcoal  to  stir  for  2-3  min  before  beginning 

addition.  With  tubes  still  In  the  ice  bath#  add  0.5  mL  of  Dextran- 
coated  charcoal  to  all  tubes  except  total  count  tubes.  Vortex  tubes 
and  Incubate  In  ice  bath  for  10  min  after  the  last  addition. 

9.  Centrifuge  all  tubes  at  3000  rpm  at  4  °C  for  10  min. 

10.  Pour  supernate  Into  scintillation  vials;  add  3.5  mL  of  scintillation 
fluid. 

11.  Count  each  vial  for  10  min  In  a  liquid  scintillation  counter.  The 
samples  are  run  with  three  pooled  (beginning#  middle#  and  end  of  run) 
rat  sera  to  assure  within-run  and  between-run  control.  Standard 
curves  are  plotted  and  values  determined  by  a  computer  program  with  a 
logit-log  plot. 


Different  volumes  of  diluted  plasma  (5#  10#  25#  and  50  pL)  were  assayed  and 
showed  good  linearity.  The  Intraassay  coefficient  of  variation  of  a  pool 
containing  10.35  pg/dL  was  5.9%  (n  =  9)  and  of  a  pool  containing  31.9  yg/dL  was 
6.6%  (n  »  9).  The  Interassay  variability  of  the  latter  pool  was  12.0%  (n  =  8). 
Accuracy  was  examined  by  adding  corticosterone#  1.25#  2.50#  5#  10#  20#  and  40 
yg/dL#  to  charcoal -treated  plasma  and  assaying  the  plasma.  The  regression  line 
was  y  =  1.009x  +  0.49#  Indicating  no  systematic  error. 

Resting  Concentrations  of  Plasma  ACTH  and  Plasma  Corticosterone.  Several 
groups  of  preliminary  experiments  were  performed  to  determine  the  real  resting 
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value  of  plasma  ACTH  and  plasma  corticosterone  In  the  cannulated  rats  (Figs.  1# 
2,  4,  5,  and  6). 

There  are  reports  that  plasma  resting  levels  of  ACTH  and  corticosterone  are 
smaller  In  decapitated  than  In  the  cannulated  rats.  However*  most  Investigators 
have  found  that  the  values  of  plasma  ACTH  or  corticosterone  obtained  from 
decapitated  rats  are  similar  to  the  values  obtained  from  cannulated  rats.  This 
finding  was  confirmed  during  preliminary  experiments  under  this  study.  Other 
reports  Indicate  that  Individually  housed*  cannulated  rats  have  slightly 
elevated  morning  levels  of  plasma  ACTH  and  corticosterone.  Evening  or  noise- 
stimulated  hormone  levels  In  these  rats  were  similar  to  those  obtained  from 
decapitated  rats.  In  this  study*  blood  sampled  between  10  AM  and  noon  from 
either  cannulated  animals  or  from  decapitated  animals  had  similar  values  of 
plasma  corticosterone. 
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III.  RESULTS  AND  ANALYSIS 


Plasma  ACTH.  Appendix  A  contains  the  data  collected  during  the  pre¬ 
radiation.  radiation,  and  postradiation  periods  for  both  the  exposure  and  sham- 
exposure  groups.  The  high  variance  displayed  by  the  data  for  the  entire 
sampling  period  Indicated  various  levels  of  physiological  activity  at  the  time 
of  blood  sampling.  Since  the  boxes  had  opaque  walls,  the  activity  of  each 
animal  prior  to  sampling  was  not  recorded.  However,  each  animal  had  sufficient 
time  (30  min  or  more)  to  return  to  basal  hormone  level  after  the  stimulation  of 
being  placed  In  the  sampling  box. 

Figures  8  and  9  present  the  raw  ACTH  concentration  data  In  scatter  diagram 
form  (the  dotted  lines  pass  through  the  mean  response  at  each  week).  Despite 
efforts  to  condition  the  animals  to  the  sampling  box  environment  prior  to 
drawing  blood  samples,  the  basal  resting  value  of  plasma  ACTH  decreased  during 
weeks  -3.  -2,  and  0.  Beyond  this  point#  the  data  displayed  a  generally  linear 
response#  factoring  out  spikes  at  weeks  7  and  10  (sham-exposure  group)  and  at 
weeks  8,  12,  and  17  (exposure  group).  The  spike  at  week  10  (sham-exposure 
group)  did  not  receive  much  weight  since  Its  value  derived  from  a  single, 
"influential,"  observation.  Unusually  high  personnel  activity  (with  attendant 
noise)  may  have  been  the  cause  of  the  remaining  spikes. 

Plasma  ACTH  concentrations  In  the  exposure  animals  did  not  appear  to  be 
greater  than  plasma  ACTH  concentrations  In  the  sham-exposure  animals  (Fig.  10). 
This  was  preliminary  evidence  Indicating  that  435-MHz  RFR  did  not  Increase  the 
resting  level  of  plasma  ACTH.  To  attach  a  numerical  probability  to  this 
conclusion,  a  statistical  analysis  was  performed  on  the  data. 

Using  multiple  linear  regression  procedures,  a  model  was  built  to  describe 
plasma  ACTH  concentration  as  a  function  of  Incident  RFR  and  time.  Terms  of  the 
polynomial  model  tested  for  the  presence  of  RFR  and  time  effects.  Various 
statistical  diagnostic  procedures.  Including  lack-of-flt  tests  and  residual 
analysis,  were  then  applied  to  the  developed  model  to  check  Its  validity. 
Appendix  B  contains  detailed  descriptions  of  the  methodology,  procedures,  and 
results  of  both  the  ACTH  and  corticosterone  data  analyses. 

The  statistical  analysis  Indicated  that  435-MHz  RFR  did  not  Increase  plasma 
concentrations  of  ACTH  In  the  exposure  animals  when  compared  to  the  sham- 
exposure  animals.  In  fact,  the  plasma  concentration  of  ACTH  In  the  exposure 
animals  remained  lower  through  most  of  the  experiment  when  compared  to  the  sham- 
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Figure  8.  ACTH  concentration  data  scatter  diagram  (sham-exposure  group). 
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Figure  9.  ACTH  concentration  data  scatter  diagram  (exposure  group). 
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Figure  10.  Mean  plasma  ACTH  concentrations  versus  time 


exposure  animals.  At  the  experiment  onset#  the  estimated  resting  value  of  ACTH 
was  99.5  pg/mL.  Further  analysis  determl'ned  that#  If  there  were  any  RFR-Induced 
effects#  they  had  to  lie  within  a  range  of  +7.37  pg/mL  (between  92.1  to  106.9) 
from  the  resting  value.  Since  resting  ACTH  concentrations  between  50  and  150 
pg/mL  are  considered  normal  In  unstressed  rats#  then  there  was  no  Indication 
that  the  long-term  RFR  exposure  produced  any  stress  as  measured  by  plasma  ACTH 
concentration. 

Plasma  corticosterone.  Appendix  G  contains  the  data  collected  during  the 
preradiation#  radiation#  and  postradiation  periods  for  both  exposure  and  sham- 
exposure  groups.  Like  ACTH#  this  data  also  displayed  a  high  variance  (for  the 
same  reasons  as  stated  previously). 

Figures  11  and  12  present  the  raw  corticosterone  data  In  scatter  diagram 
form  (the  dotted  lines  pass  through  the  mean  response  at  each  week).  The  plasma 
corticosterone  data  were  somewhat  better  behaved  than  the  ACTH  data.  Spikes 
occurred  at  weeks  7  and  12  In  the  sham-exposure  animals  (the  spike  at  week  7 
corresponds  to  that  of  ACTH  at  week  7  and  both  were  probably  due  to  the  same 
cause),  and  at  weeks  17  and  22  In  the  exposure  animals. 

Once  again,  both  plots  appeared  to  be  essentially  the  same  within 
experimental  error;  additionally#  the  RFR-exposed  animals  had  corticosterone 
concentrations  below  that  of  their  sham-exposed  counterparts  for  most  of  the 
study  (Fig.  13).  This  was  an  indication  that  435-MHz  RFR  exposure  did  not 
Increase  the  resting  value  of  plasma  corticosterone. 

The  corticosterone  statistical  analysis  Is  detailed  In  the  latter  half  of 
Appendix  B.  The  procedure  applied  to  the  data  was  Identical  to  that  applied  to 
plasma  ACTH  data  and  the  result  obtained  was  similar.  Therefore#  435-MHz  RFR 
was  not  found  to  Increase  resting  plasma  corticosterone  concentration  In  the 
exposure  group  when  compared  to  the  sham-exposure  group.  In  the  sham-exposure 
group#  plasma  corticosterone  concentrations  rose  from  a  level  of  approximately 
13.6  yg/dL  at  the  "exposure”  onset  to  approximately  17.1  ug/dL  at  the  "exposure" 
termination.  Neither  of  these  values  was  typical  of  corticosterone 
concentration  In  stressed  rats.  In  the  exposure  group#  plasma  corticosterone 
concentrations  started  at  approximately  13.6  yg/dL  at  exposure  onset#  dropped  to 
approximately  10.3  yg/dL  at  week  9  of  the  radiation  period#  and  then  rose  to 
approximately  17.9  yg/dL  at  the  exposure  conclusion.  These  values  were  also 
atypical  of  plasma  corticosterone  concentration  In  stressed  animals.  Further 
analysis  Indicated  that#  If  there  were  any  RFR-Induced  effects  on  plasma 
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Figure  12.  Corticosterone  concentration  data  scatter  diagram  (exposure  group). 
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Figure  13.  Mean  plasma  corticosterone  concentrations  versus  time. 


corticosterone  concentration,  they  had  to  lie  within  a  range  of  +  1.31  yg/dL 
from  the  resting  value.  Since  resting  corticosterone  concentrations  below  25 
yg/dL  are  considered  normal  In  unstressed  rats,  then  there  was  no  Indication 
that  the  long-term  RFR  exposure  group  had  an  elevated  mean  resting  value  of 
plasma  corticosterone.  Therefore,  435-MHz  RFR  did  not  appear  to  Induce  any 
stress  In  the  exposure  group,  as  measured  by  the  concentration  of  plasma 
corticosterone. 
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IV.  DISCUSSION 


Stress  has  been  described  as  the  "non-specific  response  of  an  organism  to 
any  demand"  [30,31].  Both  external  and  Internal  stresses  evoke  Increased  ACTH 
and  corticosterone  secretion  In  the  chronically  cannulated  rat  [323. 

Evaluation  of  the  adrenocortical  system  In  the  rat  depended  on  accurate 
determinations  of  plasma  ACTH  and  plasma  corticosterone.  Values  of  basal  levels 
of  these  hormones  varied  considerably  among  laboratories  and  depended  also  on 
the  time  when  they  were  measured.  The  differences  resulted  from  the  Increased 
sensitivity  of  new  assays  permitting  measurement  of  very  low  concentrations  of 
ACTH  and  corticosterone  that  were  not  possible  a  few  years  ago.  Conditions 
under  which  animals  were  maintained  during  blood  sampling  were  the  source  of 
other  variations  In  resting  levels  of  plasma  hormones.  Acute  elevations  of 
plasma  ACTH  and  plasma  corticosterone  due  to  handling  may  mask  experimentally 
Induced  stimulation  of  the  adrenocortical  system.  The  cannulatlon  technique 
permitted  repeated  use  of  the  same,  unanesthetized,  resting,  and  undisturbed  rat 
for  blood  sampling  and  hormone  determination,  thereby  minimizing  the  amount  of 
handling  needed  to  obtain  a  blood  sample. 

There  are  few  studies  concerned  with  Investigating  the  effects  of  long- 
lasting  stress  on  the  pituitary-adrenocortical  axis.  For  this  reason,  the  work 
of  Burchfield  et  al.  [333  was  of  particular  Interest.  Their  rats  received 
chronic  or  repeated  acute  cold  stress  during  a  period  of  up  to  3  months.  The 
authors  found  that  chronically  stressed  rats  had  elevated  resting  plasma 
corticosterone  levels,  as  much  as  seen  In  the  control  animals  during  an  acute 
stress.  On  the  other  hand,  plasma  ACTH  levels  remained  unchanged.  Thus, 
chronic  stress  lasting  3  months  Increased  the  plasma  level  of  corticosterone  In 
undisturbed,  resting  rats  from  a  value  of  8-12  ug/dL  to  23  ug/dL,  but  left  the 
resting  plasma  ACTH  level  unchanged.  The  same  experiments  demonstrated  also 
that  longer  stress  duration  led  to  higher  plasma  corticosterone  levels. 
Sakellarls  and  Vernlcos-Danell Is  [343  found  that  after  adaptation  to  chronic 
stresses,  and  despite  continued  exposure  to  the  stressor,  the  pituitary-adrenal 
system  regulates  plasma  corticosterone  concentration  at  the  prestress  levels. 
The  pituitary  secretion  returned  to  control  value  after  the  adaptation  to  stress 
was  completed. 

The  high  sensitivity  of  the  braln-pltultary-adrenal  gland  system  requires 
only  remote  sampling  of  the  blood.  The  repeated  sampling  of  blood  from  the 
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same  cannulated  rat  provided  us  with  reliable  information  regarding  resting 
patterns  of  ACTH  and  corticosterone  secretion  as  well  as  any  Increases  Induced 
by  a  long-term*  low-level  RFR  exposure.  In  other  words*  It  seemed  likely  that 
even  the  smallest  environmental  perturbation,  such  as  low-level  RFR*  would  have 
been  easily  detectable  If  It  has  any  Influence  on  the  brain-pituitary-adrenal 
axis.  Although  relocation  of  a  rat  from  the  cage  Into  the  sampling  box  30  min 
before  blood  sampling  slightly  disturbed  the  environment  of  the  rat,  such  a 
perturbation  did  not  alter  resting  plasma  ACTH  and  plasma  corticosterone  levels 
(waiting  period  and  adaptation  of  the  rat  to  the  sampling  box). 

Short-term  exposure  to  low-level  RFR  does  not  change  plasma  corticosterone 
levels  of  exposed  rats  [35*36]  although  Guy  and  associates  [37]  found  an 
elevation  of  plasma  corticosterone  the  first  time  the  blood  was  sampled  from 
RFR-exposed  rats  In  their  long-term  study.  In  the  same  study,  plasma 
corticosterone  returned  to  resting  control  levels  throughout  the  remaining  2- 
year  period.  It  was  possible  that  the  early  Increase  In  plasma  corticosterone 
of  RFR-exposed  rats  was  the  consequence  of  excitement  resulting  from  the 
animal's  early  exposure  to  the  blood  sampling  process.  Similar  early  Increases 
in  plasma  corticosterone  have  been  observed  [32]  In  which  the  animals  were  not 
fully  preconditioned  to  blood  sampling  boxes. 

As  far  as  we  could  determine,  there  Is  no  Information  available  regarding 
the  effect  of  short-term,  low-level  RFR  on  plasma  ACTH  levels.  This  Information 
would  permit  an  evaluation  of  the  pituitary-adrenal  axis  as  a  response  to  this 
environmental  condition. 

Both  plasma  ACTH  and  plasma  corticosterone  are  sensitive  Indicators  of 
various  types  of  environmental  stress  In  mammalian  systems.  In  stressful 
situations,  the  plasma  concentrations  of  both  hormones  Increase,  and  the 
increase  observed  Is  often  a  function  of  stress  Intensity  and  duration.  For 
example,  previous  studies  [33]  have  demonstrated  that  long-term  cold  stress 
Induced  In  rats  a  several-fold  Increase  In  resting  corticosterone  concentrations 
without  a  corresponding  Increase  In  ACTH  concentrations.  Johnson  et  al.  [37] 
reported  that  a  2450-MHz  pulsed-wave  RFR  environment  may  have  been  minimally 
stressful  to  rats  over  a  2-year  exposure  duration.  Their  conclusion  was  In  part 
based  on  an  elevation  of  plasma  corticosterone  levels  detected  during  their 
first  assessment  period,  followed  by  a  return  to  basal  levels  for  the  remainder 
of  the  exposure,  as  well  as  a  decrease  In  open-field  behavior  In  the  exposure 
group  during  the  first  assessment  period. 
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Our  results  Indicated  that  a  435-MHz  pulsed-wave  RFR  environment  did  not 
Increase  either  resting  ACTH  or  resting  corticosterone  concentrations  In  rats. 
Further#  a  statistical  analysis  of  the  data  Indicated  that  If  there  were  any 
RFR-induced  effects  on  resting  plasma  hormone  concentrations#  they  would  lay 
within  a  range  of  +  7.4  pg/mL  from  an  estimated  resting  concentration  of  99.5 
pg/mL  In  ACTH;  +  1.3  yg/dL  from  an  estimated  resting  concentration  of  13.6pg/dL 
In  corticosterone.  These  values  are  not  typical  of  rats  exposed  to  stress. 
Therefore#  this  study  concludes  that  a  1.0  mW/cm^  435-MHz  pulsed-wave  (1.0-ps 
pulse  width#  1-kHz  pulse  rate)  RFR  environment  did  not  Induce  any  detectable 
Increase  In  stress#  as  measured  by  resting  ACTH  and  corticosterone 
concentrations#  In  the  exposure  group  of  87  cannulated  male  Sprague-Dawley  rats 
when  compared  to  a  sham-exposure  group  of  65  cannulated  male  Sprague-Dawley 
rats. 
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APPENDIX  B 


STATISTICAL  METHODOLOGY 


The  balanced  design  of  this  experiment  (requiring  that  25  animals  from  each 
100  animal  group  be  sampled  once  every  3  weeks  for  stress  hormones)  should  have 
produced  data  easily  tested  by  balanced#  2-way  analysis  of  variance  (ANOVA) 
statistics  with  12  levels  of  factor  A  (time)  and  2  levels  of  factor  B  (RF 
radiation).  However#  data  collection  did  not  proceed  according  to  protocol  In 
that#  in  numerous  cases#  samples  were  collected  at  odd  Intervals  (Invalidating 
the  orthogonality  of  the  design)  and  the  number  of  samples  taken  per  week  varied 
above  and  below  the  25  animal  mark  (unbalancing  the  design).  These  two  factors 
combined  to  lower  the  power  of  ANOVA  statistics  (power  being  defined  as  the 
ability  to  reject  the  null  hypothesis  given  the  null  hypothesis  should  be 
rejected)  trying  to  test  the  model: 


yijk =  v  +  m  +  6j  +  T6tj + 


(B-l) 


where 


T1 

Bj 


1  jk 


=  hormone  concentration  (response), 

=  the  normal  hormone  resting  concentration# 

=  the  change  In  hormone  resting  concentration  Induced  by  RFR# 

=  the  change  In  hormone  resting  concentration  Induced  by  time, 

*  the  change  In  hormone  resting  concentration  Induced  by  the 
Interaction  between  RFR  and  time#  and 
=  noise  within  the  system  (sampling  and  assaying  errors) 


for  the  following  hypotheses: 


H0:  t0  =  Ti  =  0,  ( B-2 ) 

H^:  Tq  or  ^  0  (RFR-Induced  effects)# 

Hq:  &i  «  e2  “  •  •  •  *  Bi2  =  (B-3) 

H^:  at  least  one  Bj  f  0  (time-induced  effects)# 

Hq:  =  0#  and  (B-4) 

H^;  at  least  one  ^  0  (Interaction  between  RFR  and  time). 
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However*  examination  of  the  collected  data  suggested  an  alternative 
approach  In  that  the  data  resembled  what  might  have  been  collected  In  an 
unplanned  experiment  monitoring  over  time  the  operation  (In  this  case, 
characterized  by  resting  animal  hormone  concentrations)  of  an  established  RF 
radiation  facility.  Data  of  this  type  are  often  successfully  treated  by 
employing  linear  regression  techniques  to  develop*  build,  and  test  a  linear  (or 
Intrinsically  linear)  model  whose  parameters  can  be  used  to  predict  the  system 
response  at  various  treatment  levels.  Therefore,  we  decided  to  proceed  with  a 
regression  approach  to  data  analysis. 

Plasma  ACTH  Statistical  Analysis 

Examination  of  the  ACTH  scatter  diagrams  of  Figures  8  and  9  showed  an 
essentially  linear  response  (factoring  out  spikes  displayed  at  weeks  7,  10,  8, 
12,  and  17)  of  plasma  ACTH  versus  time.  Therefore,  the  linear  model  to  test  for 
RFR-Induced  effects  on  the  ACTH  concentration  would  Include  a  nonzero  Bg  term 
(which  provided  an  estimate  of  the  mean  hormone  resting  level  at  week  0),  and  an 
cxgz  term  to  test  for  any  biasing  of  normal  resting  levels  due  to  RFR.  Since  the 
line  traced  by  the  mean  hormone  values  In  the  scatter  diagrams  was  not 
completely  straight,  linear  terms  (3^  and  ot^z)  were  Included  In  the  .model  to 
adjust  for  a  change  In  slope  over  time,  and  quadratic  terms  ($^  and  a^z)  were 
Included  In  the  model  to  adjust  for  any  curvature  in  the  response. 

Therefore,  the  Initial  model  became: 

y  =  3g  +  3]X  +  Bnx2  +  aQz  +  aizx  +  a  nzx2,  (B-5) 

where  y  =  plasma  ACTH  concentration, 
x  *  time  (in  weeks),  and 

z  =  a  categorical  variable  with  value  0  for  animals  In  the  sham- 
exposure  group  and  value  1  for  animals  In  the  exposure  group. 

Raw  data  from  the  ACTH  spreadsheet  (Included  In  Appendix  A)  were  put  on 
computer  file,  and  a  Statistical  Analysis  System  (SAS)  formatting  program 
(Included  In  Appendix  C)  was  prepared  to  read  the  data  file  and  perform  the 
desired  statistical  tests  on  the  model. 

The  first  test  Identified  terms  within  the  general  model  which  contributed 
the  least  toward  forming  a  statistically  significant  regression.  These 
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procedures  were  used  In  combination  with  an  Initial  regression  on  the  general 
model  (not  Included)  to  evaluate  the  statistical  significance  of  terms  modelling 
the  ACTH  concentration  time  dependency  and  terms  modelling  the  RFR-Induced 
effects  on  ACTH  concentration.  Two  types  of  model  "building”  procedures  were 
employed:  forward  stepwise  regression  and  maximum  R2  regression.  Forward 

stepwise  regression  produced  a  "best"  model  (In  that  It  minimized  the  number  of 
model  predictor  variables)  In  which  all  the  Included  terms  were  statistically 

O 

significant  at  a  level  of  0.15.  Maximum  R^  regression  approached  model  building 
from  a  somewhat  different  perspective;  It  first  built  the  "best"  one-predictor 
model  (in  that  It  maximized  R2,  the  percentage  of  variation  within  the  data 
explained  by  the  regression),  then  built  the  "best"  two-predictor  model,  the 
"best"  three-predictor  model,  etc.,  until  it  arrived  at  the  general  model.  All 
these  tests  Indicated  that  the  model  which  best  fit  the  data  was 

y  =  Bq  +  31x  +  B^x2  +  “jxz,  ( B— 6 ) 

where 

BQ  =  99.50, 

=  1.79, 

Bn  =  0.12, 

a  -L  =  -0.68, 

and  the  other  variables  were  as  defined  previously.  The  outputs  of  both 
procedures  are  Included  In  Appendix  D.  The  absence  of ag  Indicated  that  RFR  did 
not  produce  a  detectable  effect  on  the  Intercept  of  the  model,  and  therefore  did 
not  bias  the  ACTH  concentration  of  the  exposure  group  when  compared  to  the  sham- 
exposure  group.  Also,  RFR  did  not  produce  a  detectable  effect  on  the  curvature 
of  the  exposure  group's  ACTH  concentration  when  compared  to  the  sham-exposure 
group. 

The  sham-exposure  group  did  display  a  difference,  when  compared  to  the 
exposure  group.  In  the  overall  time  response  of  hormones.  In  both  groups,  ACTH 
concentrations  at  the  onset  of  exposure  were  somewhat  high  (106  pg/mL  for  sham- 
exposure  animals,  108.1  pg/mL  for  exposure  animals),  dropped  Into  the  middle  of 
the  exposure  duration  (93.1  pg/mL  at  week  11  for  sham-exposure  animals,  87.1 
pg/mL  at  week  14  for  exposure  animals),  and  then  rose  at  exposure  termination 
(151.6  pg/mL  5  weeks  after  exposure  termination  for  sham-exposure  animals,  131.8 
pg/mL  5  weeks  after  exposure  termination  for  exposure  animals).  The  differences 


42 


In  predicted  ACTH  concentration  at  the  end  of  the  exposure  were  enlightening. 
Basically,  the  trend  (modelled  by  the  negative  a  ^xz  term)  was  that  ACTH 
concentrations  In  the  exposure  group  remained  lower  than  ACTH  concentrations  In 
the  sham-exposure  group  for  the  exposure  duration.  Both  ranges  (93.1  to  151.6 
pg/mL  In  sham-exposure  animals,  87.1  to  131.8  pg/mL  In  exposure  animals)  were 
still  well  within  the  normal  range  of  plasma  ACTH  In  non-stressed  male  Sprague- 
Dawley  rats  (approximately  75-175  pg/mL).  Therefore,  these  results  Indicated 
that  435-MHz  RFR  did  not  Induce  an  elevation  of  ACTH  concentration  In  the 
exposure  group. 

The  just  mentioned  conclusions  could  be  accepted  only  after  the  assumptions 
underlying  the  model  used  to  draw  the  conclusions  were  verified.  These 
assumptions  Included  no  model  lack-of-flt,  NID(0»  a  ^)  residual  distribution,  and 
no  model  multlcoll Inearlty. 

Since  multiple  observations  of  ACTH  concentration  were  taker  for  the  weeks 
containing  data.  It  was  possible  to  perform  a  model  lack-of-flt  test  on  the 
regression.  The  lack-of-flt  test  Involved  breaking  the  sum-of-squares  error 
from  the  regression  Into  two  components:  sum-of-squares  pure  error, 
representing  the  actual  variation  due  to  the  sampling  and  assaying  process;  and 
sum-of-squares  lack-of-flt,  representing  the  variation  due  to  the  difference 
between  the  mean  value  at  one  week  when  compared  to  the  fitted  value  at  the  same 
week.  A  test  statistic  was  then  computed  comparing  the  sum-of-squares  lack-of- 
flt  to  the  sum-of-squares  pure  error;  sufficiently  high  values  of  the  test 
statistic  Indicated  model  lack-of-flt. 

Sum-of-squares  error  was  obtained  from  the  ANOVA  table  produced  In  the 
regression  procedure  output.  Sum-of-squares  pure  error  was  obtained  by 
analyzing  the  experiment  from  2-way,  fixed  effects  ANOVA  viewpoint.  The  sum-of- 
squares  lack-of-flt  was  then  computed  from  the  difference  of  sum-of-squares 
error  minus  sum-of-squares  pure  error.  Calculations  to  compute  the  critical 
value  Fq  are  detailed  In  Appendix  E. 

The  computed  test  statistic  Fq  exceeded  the  critical  value,  thereby 
Indicating  significant  lack-of-flt.  Normally,  this  result  would  be  somewhat 
disturbing  since  It  would  require  refitting  the  model  using  transformed  rather 
than  raw  data  values.  In  fact,  transformation  of  the  dependent  variable  y  was 
definitely  undesirable,  since  the  residual  plots  Indicated  that  the  residuals  of 
y  (using  the  revised  model)  conformed  to  the  NID(0,  o^)  requirement. 
Additionally,  the  transformation  of  the  predictor  variables  x  and  x^  to  yield  a 
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model  displaying  no  lack-of-flt,  although  theoretically  possible,  would  be  a 
long  and  time-consuming  process. 

Fortunately,  the  experimental  design  helped  compensate  for  the  model  lack- 
of-flt  deficiency.  First  of  all.  the  lack-of-flt  was  comparatively  small. 
Under  optimal  conditions  (lack-of-flt  statistically  Insignificant),  both  the 
mean  square  error  and  the  mean  square  pure  error  estimate  the  population 
variance.  If  there  Is  a  lack-of-flt.  the  mean  square  pure  error  estimates  the 
population  variance,  while  the  mean  square  error  estimates  the  population 
variance  plus  a  bias  term.  From  the  ANOVA  (regression  and  GLM)  tables,  the 
tabulated  values  for  MSg  and  MSpe  were  1247.19  and  1132.31  respectively.  These 
values  corresponded  to  sample  standard  deviations  of  35.32  pg/mL  and  33.65 
pg/mL.  Thus,  although  the  lack-of-flt  was  statistically  significant.  It  was 
also  practically  Insignificant.  In  other  words,  the  development  of  an 
alternative  model  displaying  no  lack-of-flt  would  yield  essentially  (within  1  or 
2%)  the  same  results  as  the  present  model  displaying  lack-of-flt.  Rather  than 
identify  an  alternate  model  (which  would  not  be  that  much  better  a  predictive 
tool  than  the  model  currently  being  used),  we  decided  to  proceed  with  the 
stepwise  model  and  modify  the  significance  of  the  tests  to  compensate  for  model 
lack-of-flt.  Thus,  all  a  values  are  somewhat  higher  than  they  should  be.  and 
the  confidence  intervals  established  are  somewhat  wider  than  Indicated  In  the 
appendix  tables. 

The  next  step  In  determining  model  accuracy  Involved  examining  the  residual 
and  partial  residual  plots  to  verify  the  least-squares  regression  assumption 
that  the  model  errors  were  NID(O.o^),  Confirming  this  assumption  defended  the 
use  of  F  tests  to  determine  the  statistical  significance  of  the  parameters,  and 
validated  the  statistics  which  produced  the  tables  listing  confidence  Intervals 
of  the  ACTH  concentrations.  A  number  of  residual  plots  suggested  themselves 
Immediately:  residuals  versus  time,  residuals  versus  predicted  value  of  ACTH 

concentration,  residuals  versus  animal  case  number,  studentlzed  residuals  versus 
the  above  three,  and  partial  residual  plots  corrected  for  the  model  terms  £q» 
B^»  B^.  and  a  j_.  Since  the  order  of  each  sample's  assay  was  unavailable,  there 
was  no  way  to  test  for  correlation  using  either  the  Durbln-Watson  or  runs  test. 

Examination  of  the  residual  plots  yielded  no  discernible  pattern  In  the 
distribution  of  the  residuals,  thus  confirming  the  NID(O.o^)  hypothesis. 
However,  four  outliers  at 
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animal  #06 
animal  #40 
animal  #44 
animal  #47 


week  18  [ACTH]  =  290  pg/mL 
week  16  [ACTH]  =  265  pg/mL 
week  7  [ACTH]  =  287  pg/mL 
week  7  [ACTH]  =  299  pg/mL 


were  detected  using  the  criterion  of  Cook’s  distance  (each  of  the  above 
observations  had  a  Cook’s  D  of  approximately  1%,  making  their  contribution  to 
the  model  Inappropriately  large  considering  the  size  of  the  data  set).  These 
outliers  were  most  likely  the  result  of  assaying  or  reporting  error*  and  were 
considered  sufficiently  anomalous  to  be  rejected  from  the  data  set.  Residual 
plots  were  then  regenerated  (see  Appendix  F)  and  checked  for  patterns.  The  new 
residual  plots  In  turn  did  not  challenge  the  ~NID(0,o2)  assumption  Inherent  In 
least-squares  model  fitting. 

To  complete  the  analysis*  diagnostics  to  check  for  model  multlcoll Inearlty 
and  correlation  between  the  terms  were  employed.  Examination  of  the  listed 
condition  numbers  and  matrix  eigenvalues  (being  provided  under  separate  cover) 
detected  no  troublesome  values.  This  analysis  Indicated  that  the  model  did  not 
display  a  significant  degree  of  multlcoll Inearlty.  Similarly*  examination  of 
the  correlation  matrix  showed  that  correlation  between  the  estimated  values  ofS 
were  all  within  tolerable  limits.  The  highest  degree  of  correlation  was  between 
the  x  and  the  x2  term,  which  often  occurs  when  using  a  polynomial  model  In 
linear  regression. 

For  future  reference,  and  for  the  sake  of  completeness*  tables  listing 
animal  case  number*  observations  (If  taken)  at  each  week*  predicted  value  of 
ACTH  concentration,  standardized  error  of  prediction,  95%  confidence  Intervals 
on  the  mean  value  of  the  ACTH  concentration*  and  residuals  were  prepared,  as 
were  tables  containing  animal  case  number*  regular  and  studentlzed  residual 
values*  a  graphical  display  of  student  residual  values*  and  Influence  statistics 
(such  as  Cook's  D).  These  tables  were  used  to  detect  both  outliers  and 
Influential  data  points  In  the  ACTH  data  set. 

To  arrive  at  a  conservative  estimate  of  the  minimum  change  due  to  RFR  In 
resting  ACTH  concentrations  which  this  protocol  was  capable  of  detecting*  the 
value  of  the  operating  curve  parameter  $  g  corresponding  to  the  time  factor  (B) 
discussed  at  the  beginning  of  the  statistical  methodology  was  calculated.  This 
parameter  was  given  by 


$ 


2 

B 


naD2 
Tto1  ’ 


(B-7) 
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where  n  =  number  of  replications  per  cell  =  40» 
a  =  number  of  levels  of  factor  A  =  12# 
b  =  number  of  levels  of  factor  B  =  2# 
o -  population  variance#  and 
D2  =  detection  threshold. 

Substituting  In  values  for  a#  b»  n»  and  the  MSp0  as  an  estimate  of  a  2 
provided  an  operating  curve  parameter  of 

*  B  =  0.3255  D  (B-8) 

To  obtain  a  value  of  *  from  the  operating  curve#  the  type  I  risk  a  and  type  II 
risk  6  were  set  to  0.05  and  0.10#  respectively.  Then#  the  value  of  $  was  read 
from  the  fixed  effects  ANOVA  curve  with  v  ^  =  i  and  v  2  ■  936.  This  value  was 

*B  =  2.4  (B-9) 

Note  that  the  numerator  degrees  of  freedom#  v and  the  denominator  degrees 
of  freedom  ^2#  were  calculated  from  the  equations 

v  2  =  b-1#  and  (B-10) 

v2  =  ab(n-l) .  (B-ll) 

The  detection  level  was  therefore 

Db  =  7.37  pg/mL.  (B-12) 

Thus#  this  protocol  conservatively  was  able  to  detect  an  Increase  In 
resting  plasma  ACTH  concentrations  of  7.37  pg/mL  approximately  90%  of  the  time. 

Ej  isma  Corticosterone  Statistical  Analysis 

Examination  of  the  corticosterone  scatter  diagrams  of  Figures  11  and  12 
showed  a  situation  similar  to  that  of  ACTH.  The  mean  corticosterone  values# 
excepting  spikes  at  weeks  7#  12#  17  and  22#  followed  an  essentially  l*iear 
response  over  time.  Therefore#  the  linear  model 

y  “  ®o  +  ^lx  +  +0tOz  +  °12X  +  a  nzx^'  (B— 13 ) 
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where 


y  =  plasma  corticosterone  concentration* 
x  =  time  (In  weeks)*  and 

z  =  a  categorical  variable  set  to  0  for  the  sham-exposure  animals  and 
1  for  the  exposure  animals* 

was  again  employed.  Data  from  the  corticosterone  spreadsheet  (Included  In 
Appendix  G)  were  put  Into  a  new  file  and  a  second  SAS  formatting  program 
(Included  In  Appendix  H)  was  prepared  to  analyze  the  data. 

The  model  Indicated  by  the  forward  and  maximum  stepwise  regression 
procedures  was 

y  =  Bq  +  B]_x  +  iijxz  +  (B-14) 

where  Bg  =  13.55* 

Bj.  =  0.15, 

=  -0.83,  and 
a  u  =  0.036, 

with  the  x,  y»  and  z  variables  as  defined  previously.  The  outputs  of  both 
procedures  are  Included  In  Appendix  I.  The  absence  of  “g  Indicated  that  RFR  did 
not  Induce  a  detectable  change  In  the  Intercept  of  the  model  and  therefore  did 
not  bias  the  corticosterone  concentration  of  the  exposure  group  when  compared  to 
the  sham-exposure  group. 

In  the  case  of  the  sham-exposure  group*  the  resting  plasma  corticosterone 
concentration  slowly  rose  during  the  exposure  from  a  level  of  13.55  ug/dL  at  the 
study  onset  to  17.08  ^g/dL  at  week  24  of  the  exposure  schedule.  Both  of  these 
values  were  well  within  the  normal  resting  range  of  plasma  corticosterone  In 
unstressed  male  Sprague-Dawley  rats.  The  exposure  group  displayed  a  small 
amount  of  curvature  when  compared  to  the  sham-exposure  group,  with 
corticosterone  concentrations  of  13.55  ^g/dL  at  exposure  onset*  dropping  to 
10.32  ^g/dL  at  week  9,  and  then  rising  to  17.93  ^g/dL  at  week  24  of  the  exposure 
schedule.  Once  again,  these  values  were  well  within  the  normal  resting  range  of 
plasma  corticosterone  In  unstressed  male  Sprague-Dawley  rats.  Therefore*  from 
the  above  model*  RFR  did  not  Induce  an  elevation  of  plasma  corticosterone  In  the 
exposure  group.  The  trend  displayed  In  the  data  was  that  sham-exposure 
corticosterone  concentrations  gradually  Increased  over  time*  the  Increase 
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remaining  well  within  the  boundaries  of  plasma  corticosterone  In  unstressed 
rats.  The  exposure  group's  corticosterone  levels  tended  to  fall  at  the  exposure 
C!  ^et,  and  consistently  remained  lower  than  those  of  the  sham-exposure  group 
until  exposure  termination*  when  the  two  group's  concentrations  were  essentially 
equal.  Although  the  two  groups  did  not  behave  In  an  identical  manner*  there  was 
no  statistical  basis  to  Indicate  that  RFR  Induced  a  stress  which  manifested 
Itself  as  a  difference  In  plasma  corticosterone  concentration  between  the 
exposure  and  sham-exposure  groups. 

It  was  then  necessary  to  check  the  model  lack-of-fit*  the  residual 
distribution*  and  the  amount  of  multlcoll inearlty  present  in  the  model.  First* 
the  model  was  checked  for  lack-of-fit.  Sum-of-squares  error  from  the  regression 
was  92991.99  with  1297  degrees  of  freedom  and  sum-of-squares  pure  error  was 
81974.68  with  1254  degrees  of  freedom.  Therefore*  sum-of-squares  lack-of-fit 
was  11017.31  with  43  degrees  of  freedom*  and  Fg  was  calculated  to  be  3.9195  (see 
Appendix  J  for  details  of  this  calculation). 

As  In  ACTH,  the  model  derived  for  corticosterone  displayed  a  significant 
lack-of-fit.  Transformation  of  the  response  variable  y  was  undesirable  since 
the  residuals  of  y  (generated  with  the  given  model  predictors)  conformed  to  the 
NID(0,a2)  hypothesis.  Transformation  of  the  x  values*  while  possible*  would 
remain  a  rather  long  and  time-consuming  process. 

This  problem  was  solved  the  same  way  as  with  ACTH.  Since  the  lack-of-fit 
was  comparatively  small  (with  MSg  =  71.70  and  MSpe  =  65.37  as  estimates  of  the 
corticosterone  variance)*  the  difference  In  predicted  responses  between  a  model 
displaying  Insignificant  lack-of-fit  and  this  model  (displaying  significant  but 
small  lack-of-fit)  would  only  amount  to  values  differing  from  one  another  by  1 
to  2%.  Therefore,  an  analysis  using  the  model  that  displayed  some  lack-of-fit 
was  used.  This  lack-of-fit  was  then  compensated  for  by  altering  the 
significance  levels  a  for  the  conclusions  drawn  from  the  analysis. 

Appendix  K  contains  the  residual  plots  generated  from  the  data  set  once  the 
outliers  were  removed.  The  plots  all  showed  no  discernible  pattern*  and  thereby 
did  not  void  the  assumption  that  the  model  e^s  were  distributed  as  normal  random 
variables.  Tests  for  Independence  were  not  possible  since  the  order  of  each 
sample's  assay  was  unavailable.  Therefore*  from  the  given  residuals  and  data, 
there  was  no  reason  to  reject  the  assumption  that  the  residuals  were  distributed 
NID(0»o2) . 
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The  outliers  detected  by  Cook's  distance  and  removed  from  the  data  set  were 
as  follows: 


animal 

#  14 

week 

14 

[cortlco] 

= 

38 

pg/dL 

animal 

#  23 

week 

4 

[cortlco] 

= 

48 

pg/dL 

animal 

#  54 

week 

2 

[cortlco] 

= 

49 

pg/dL 

animal 

#  55 

week 

2 

[cortlco] 

= 

57 

pg/dL 

animal 

#124 

week 

2 

[cortlco] 

= 

40 

pg/dL 

animal 

#125 

week 

1 

[cortlco] 

= 

41 

pg/dL 

animal 

#143 

week 

-3 

[cortlco] 

= 

38 

pg/dL 

animal 

#178 

week 

10 

[cortlco] 

= 

44 

pg/dL 

As  was  the  case  with  ACTH.  these  outliers  were  most  likely  the  result  of 
assaying  or  recording  errors,  and  were  considered  sufficiently  anomalous  to  be 
rejected  from  the  data  set. 

Examination  of  model  multlcol linearity  diagnostics  (provided  under  separate 
cover)  did  not  indicate  any  significant  model  multlcoll Inearlty.  The  values  for 
the  condition  Indices  were  all  within  the  [1.10]  range  considered  tolerable. 
The  highest  correlation  In  the  model  (-0.7784)  was  between  the  parameters  ot  y 
and  This  was  not  unexpected,  since  linear  and  quadratic  terms  often 
display  a  degree  of  correlation  In  polynomial  regression  situations. 
Additionally,  this  correlation  was  well  within  the  +  0.95  correlation  criteria 
which  often  Indicates  multlcoll Inearlty  problems.  Therefore,  we  concluded  that 
multlcoll  Inearlty  and  correlation  did  not  pose  serious  challenges  to  the 
validity  of  the  obtained  model. 

For  future  reference,  and  for  the  sake  of  completeness,  tables  listing 
animal  case  number,  observations  (If  taken)  at  each  week,  predicted  value  of 
corticosterone  concentration,  standardized  error  of  prediction.  95%  confidence 
Intervals  on  the  mean  value  of  the  corticosterone  concentration,  and  residuals 
were  prepared,  as  were  tables  containing  animal  case  number,  regular  and 
studentlzed  residual  values,  a  graphical  display  of  student  residual  values,  and 
Influence  statistics  (such  as  Cook's  D).  These  tables  were  used  to  detect  both 
outliers  and  Influential  data  points  In  the  corticosterone  data  set. 

To  arrive  at  a  conservative  estimate  of  the  minimum  change  due  to  RFR  In 
resting  plasma  corticosterone  concentrations  which  this  protocol  was  capable  of 
detecting,  a  calculation  was  made  for  the  value  of  the  operating  curve  parameter 
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$g  corresponding  to  the  RFR  factor  discussed  at  the  beginning  of  this  Section. 
Using  the  equation 


-2  _  naD2 

"  2EE2’  (B-15) 


where 


n  =  number  of  replications  per  cell  =  40# 
a  =  number  of  levels  of  factor  A  =  12# 
b  =  number  of  levels  of  factor  B  =  2# 
o2  =  population  variance#  and 
=  detection  threshold 


and  the  MS  as  an  estimate  of  ®2  yielded  a  value  of 

pe 


$  B  =  1.S357 


(B-16) 


To  obtain  a  value  of  $  from  the  ANOVA  operating  curve,  the  type  I  risk  a  was 
set  to  0.05#  the  type  II  risk  6  was  set  to  0.10#  and  the  degrees  of  freedom  In 
the  denominator  was  set  to  the  error  degrees  of  freedom.  Then,  reading  from  the  i 
operating  curve#  It  was  observed  that 


$  =  2.4. 


(B-17) 


Therefore#  a  conservative  estimate  was  that  the  protocol  could  detect  an  1.31 
ug/dL  Increase  In  resting  plasma  corticosterone  concentrations  approximately  90% 
of  the  time. 

We  thank  Dr.  Russell  G.  Helkes  of  Georgia  Tech's  Department  of  Industrial 
and  Systems  Engineering  for  his  assistance  In  developing  the  statistical 
methodology  of  this  appendix.  We  also  thank  Sandra  L.  Barnes  for  her  assistance 
In  analyzing  the  corticosterone  data  set. 
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APPENDIX  C 


ACTH  SAS  FORMATTING  PROGRAM 
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1  SAS(R)  LOG  CMS  SAS  5.16  VM/CMS  CMS  USER  QSECLSB 

NOTE:  COPYRIGHT  (C)  1984,1986  SAS  INSTITUTE  INC.,  CARY,  N.C.  27511,  U.S.A. 
NOTE:  CMS  SAS  RELEASE  5.16  AT  GEORGIA  INSTITUTE  OF  TECHNOLOGY  (03559001). 

NOTE:  CPUID  VERSION  -  FF  SERIAL  -  012242  MODEL  -  4381  . 

NOTE:  SAS  OPTIONS  SPECIFIED  ARE: 

LEAVE=0 

1  DATA  TESTA; 

2  CMS  FILEDEF  X  DISK  ACTH  DAT  A; 
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Al 

( 

BLKSIZE 

141 

RECFM 

VBA 

LRECL 

133 

11  ARRAY  WEEK  {33}  WKN3  WKN2  MISSN1  WK0-WK24  MISS25  WKP2  MISS27  MISS28  WKP5 

12  KEEP  X  XSQR  Y  2  XZ  XSQRZ  CASE; 

13  INFILE  X; 

14  INPUT  CASE  1-3 


15 

WKN3 

5-7 

16 

WKN2 

9-11 

17 

WKO 

13-15 

18 

WK1 

17-19 

19 

WK2 

21-23 

20 

WK3 

25-27 

21 

WK4 

29-31 

22 

WK5 

33-35 

23 

WK6 

37-39 

24 

WK7 

41-43 

25 

WK8 

45-47 

26 

WK9 

49-51 

27 

WK10 

53-55 

28 

WK11 

57-59 

29 

WK12 

61-63 

30 

WK13 

65-67 

31 

WK14 

69-71 

32 

WK15 

73-75 

33 

WK16 

77-79 

34 

WK17 

81-83 

35 

WK18 

85-87 

36 

WK19 

89-91 

37 

WK20 

93-95 

38 

WK21 

97-99 

39 

WK22 

101-103 

40 

WK23 

105-107 

41 

WK24 

109-111 

42 

WKP2 

113-115 

43 

WKP5 

117-119 

44 

t 

45 

MISSN1-. ; 

46 

MISS25-. ; 

47 

MISS27-. ; 

48 

MISS28-. ; 

49  IF  CASE  <  100  THEN  Z  -  0; 

50  IF  CASE  >«  100  THEN  Z  -  1; 


52 


2 


SAS(R)  LOG  CMS  SAS  5.16 


VM/ CMS  CMS  USER  QSECLSB 


51  IF  Z-l  THEN  CASE-CASE- 100; 

52  DO  I  -  1  TO  33; 

53  X  -  1-4;  XSQR  -  X*X;  XZ  -  X*Z;  XSQRZ  -  X*X*Z;  Y  -  WEEK  {l} ; OUTPUT; 

54  END; 


NOTE:  INFILE  X  IS  FILE  ACTH  DAT  Al 
NOTE:  151  LINES  WERE  READ  FROM  INFILE  X. 

NOTE:  DATA  SET  WORK. TESTA  HAS  4983  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.82  SECONDS  AND  248K. 

55  PROC  CONTENTS; 

NOTE:  THE  PROCEDURE  CONTENTS  USED  0.22  SECONDS  AND  376K  AND  PRINTED  PAGES  1  TO  2. 

56  PROC  PRINTTO  NEW  UNIT-20; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 

57  PROC  SORT  OUT-SCTR; 

58  BY  Z  X  Y; 

NOTE:  DATA  SET  WORK.SCTR  HAS  4983  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  PROCEDURE  SORT  USED  1.11  SECONDS  AND  6968K. 

59  PROC  SUMMARY; 

60  BY  Z  X; 

61  VAR  Y; 

62  OUTPUT  OUT-OVLMN  MEAN-MEAN; 

NOTE:  THE  DATA  SET  WORK.OVLMN  HAS  66  OBSERVATIONS  AND  5  VARIABLES. 

NOTE:  THE  PROCEDURE  SUMMARY  USED  0.85  SECONDS  AND  440K. 

63  DATA  SACTH; 

64  SET  SCTR  OVLMN; 

65  BY  Z; 

NOTE:  DATA  SET  WORK. SACTH  HAS  5049  OBSERVATIONS  AND  10  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.83  SECONDS  AND  312K. 

66  PROC  PLOT  NOLEGEND  DATA-SACTH; 

67  BY  Z; 

68  PLOT  MEAN*X- 1 X 1  Y*X-'.'  /  VAXIS-0  TO  225  BY  25  OVERLAY; 

69  TITLE  'ACTH  SCATTER  DIAGRAM'; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.58  SECONDS  AND  440K  AND  PRINTED  PAGES  3  TO  4. 

70  PROC  PRINTTO  NEW  UNIT-21; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 

71  PROC  PLOT  NOLEGEND  DATA-SACTH; 

72  PLOT  MEAN*X- ' X '  /  VAXIS-0  TO  225  BY  25; 

73  TITLE  'Mean  Plasma  ACTH  Concentration  Versus  Time'; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.23  SECONDS  AND  440K  AND  PRINTED  PAGE  5. 

74  PROC  PRINTTO  NEW  UNIT-22; 

75  TITLE  'ACTH  ANALYSIS'; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 
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3  SAS(R)  LOG  CMS  SAS  5.16  VM/CMS  CMS  USER  QSECLSB 

76  PROC  DATASETS; 


77 

•  —  j 

LIST  OF 

MEMBERS 

BEFORE  UPDATE  OF  DIRECTORY. 

NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

OVLMN 

/DATA 

66 

1 

SACTH 

/DATA 

5049 

1 

SCTR 

/DATA 

4983 

1 

TESTA 

/DATA 

4983 

1 

77 

DELETE 

SCTR; 

78 

DELETE 

OVLMN; 

79 

DELETE 

TESTA; 

LIST  OF 

MEMBERS 

AFTER  UPDATE  OF 

DIRECTORY. 

NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

SACTH 

/DATA 

5049 

1 

NOTE:  THE  PROCEDURE  DATASETS  USED  0.13  SECONDS  AND  440K. 

80  PROC  STEPWISE  DATA-SACTH; 

81  MODEL  Y  -  X  XSQR  Z  XZ  XSQRZ  /  STEPWISE  MAXR; 

NOTE:  THE  PROCEDURE  STEPWISE  USED  0.82  SECONDS  AND  440K  AND  PRINTED  PAGES  6  TO  9. 

82  PROC  PRINTTO  NEW  UNIT-23; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 

83  PROC  REG; 

84  MODEL  Y  -  X  XSQR  XZ  /  PARTIAL; 

85  ID  CASE; 

NOTE:  ACOV  AND  SPEC  OPTION  ONLY  VALID  WITH  RAWDATA 

NOTE:  THE  PROCEDURE  REG  USED  2.83  SECONDS  AND  632K  AND  PRINTED  PAGES  10  TO  14. 

86  PROC  PRINTTO  NEW  UNIT-24; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 

87  PROC  GLM; 

88  CLASS  X  Z; 

89  MODEL  Y  -  X  X*X  X*Z; 

NOTE:  THE  PROCEDURE  GLM  USED  5.17  SECONDS  AND  1016K  AND  PRINTED  PAGES  15  TO  16. 

90  PROC  PRINTTO  NEW  UNIT-25; 

91  * - 

92  * 

93  *  to  obtain  tables  listing  the  variance  inflation  factors, 

94  *  influence  statistics,  and  tolerances,  the  following  SAS 

95  *  statements  were  used  in  this  partition: 

96  * 

97  *  PROC  REG; 

98  *  MODEL  Y  -  X  XSQR  XZ  /  TOL  VIF  INFLUENCE; 

99  *  ID  CASE; 

100  *  OUTPUT  OUT-RACTH  P-PREDICT  R-RESID  STUDENT-STUDENT; 

101  * 

102  * - 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.04  SECONDS  AND  312K. 

103  PROC  REG; 

104  MODEL  Y  -  X  XSQR  XZ  /  I  SSI  SS2  STB  COVB  CORRB  SEQB  COLLIN 

105  COLLINOINT  ACOV  P  R  CLM; 


* 

* 

★ 

★ 

* 

* 

•k 

* 

* 

* 

• 
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4  SAS(R)  LOG  CMS  SAS  5.16  VM/CMS  CMS  USER  QSECLSB 

106  ID  CASE; 

107  OUTPUT  OUT-RACTH  P-PREDICT  R-RESID  STUDENT-STUDENT; 

NOTE:  THE  DATA  SET  WORK . RACTH  HAS  5049  OBSERVATIONS  AND  13  VARIABLES. 

NOTE:  THE  PROCEDURE  REG  USED  11.26  SECONDS  AND  632K  AND  PRINTED  PAGES  17  TO  123. 

108  PROC  PRINTTO  NEW  UNIT-26; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 

109  PROC  PLOT  DATA-RACTH; 

110  PLOT  RES ID*X-'*' ; 

111  PLOT  RESID*PREDICT- ' * '  /  HAXIS  -  80  TO  140  BY  5; 

112  PLOT  STUDENT*X- ' * ' ; 

113  PLOT  STUDENT*PREDICT» ' * '  /  HAXIS  -  80  TO  140  BY  5; 

114  TITLE  ' ACTH  RESIDUAL  PLOTS'; 

NOTE:  THE  PROCEDURE  PLOT  USED  2.06  SECONDS  AND  440K  AND  PRINTED  PAGES  124  TO  127 

115  PROC  PRINTTO  NEW  UNIT-27; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  312K. 

116  PROC  PLOT  DATA-RACTH; 

117  BY  Z; 

118  PLOT  RESID*CASE- ' * '  /  HAXIS  -  0  TO  90  BY  5; 

119  PLOT  STUDENT*CASE- ' * '  /  HAXIS  -  0  TO  90  BY  5; 

120  TITLE  'ACTH  RESIDUAL  PLOTS'; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.67  SECONDS  AND  440K  AND  PRINTED  PAGES  128  TO  131 
NOTE:  SAS  USED  6968K  MEMORY. 

NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 

PO  BOX  8000 

CARY,  N.C.  27511-8000 


APPENDIX  D 

STEPWISE  AND  MAXIMUM  R2  REGRESSION  PROCEDURES 
USED  TO  BUILD  ACTH  MODEL 
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SUMMARY  OF  STEPWISE  REGRESSION  PROCEDURE  FOR  DEPENDENT  VARIABLE 
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BOUNOS  ON  CONOITION  NUMBER:  8.573899.  62.65483 


THE  ABOVE  MODEL  IS  THE  BEST  3  VARIABLE  MODEL  FOUNO . 
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APPENDIX  E 


ACTH  LACK-OF-FIT  TEST 
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GENERAL  LINEAR  MODELS  PROCEDURE 


this  term  contains  both  sum-of-squares  pure  error 
and  sum-of-squares  lack-of-fit. 


APPENDIX  F 


ACTH  RESIDUAL  PLOTS 
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ACTH  RESIDUAL  PLOTS 
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ACTH  RESIDUAL  PlUIS 
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RAN  CORTICOSTERONE  DATA  SPREADSHEETS 
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CORTICOSTERONE  SAS  FORMATTING  PROGRAM 
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1 


SAS(R)  LOG  CMS  SAS  5.16 


VM/ CMS  CMS  USER  QSECLSB 


NOTE:  COPYRIGHT  (C)  1984,1986  SAS  INSTITUTE  INC.,  CARY,  N.C.  27511,  U.S.A. 
NOTE:  CMS  SAS  RELEASE  5.16  AT  GEORGIA  INSTITUTE  OF  TECHNOLOGY  (03559001). 

NOTE:  CPUID  VERSION  -  FF  SERIAL  -  012242  MODEL  -  4381  . 

NOTE:  SAS  OPTIONS  SPECIFIED  ARE: 

LEAVE-0 

1  DATA  TESTC; 

2  CMS  FILEDEF  X  DISK  CORTICO  DAT  A; 

3  CMS  FILEDEF  20  DISK  CORTICOO  LISTING  A1  (  BLKSI2E  141  RECFM  VBA  LRECL  133; 

4  CMS  FILEDEF  21  DISK  C0RTIC01  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

5  CMS  FILEDEF  22  DISK  CORTIC02  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

6  CMS  FILEDEF  23  DISK  C0RTIC03  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

7  CMS  FILEDEF  24  DISK  C0RTIC04  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

8  CMS  FILEDEF  25  DISK  C0RTIC05  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

9  CMS  FILEDEF  26  DISK  C0RTIC06  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

10  CMS  FILEDEF  27  DISK  C0RTIC07  LISTING  Al  (  BLKSIZE  141  RECFM  VBA  LRECL  133; 

11  ARRAY  WEEK  {33}  WKN3  WKN2  MISSN1  WK0-WK24  MISS25  WKP2  MISS27  MISS28  WKP5 ; 

12  KEEP  X  XSQR  Y  Z  XZ  XSQRZ  CASE; 

13  INFILE  X*. 

14  INPUT  CASE  1-3 

15  WKN3  5-6 

16  WKN2  8-9 

17  WKO  11-12 

18  WKl  14-15 

19  WK2  17-18 

20  WK3  20-21 

21  WK4  23-24 

22  WK5  26-27 

23  WK6  29-30 

24  WK7  32-33 

25  WK8  35-36 

26  WK9  38-39 

27  WK10  41-42 

28  WKl  1  44-45 

29  WK12  47-48 

30  WK13  50-51 

31  WK14  53-54 

32  WK15  56-57 

33  WK16  59-60 

34  WKl 7  62-63 

35  WK18  65-66 

36  WK19  68-69 

37  WK20  71-72 

38  WK21  74-75 

39  WK22  77-78 

40  WK23  80-81 

41  WK24  83-84 

42  WKP2  86-87 

43  WKP5  89-90 

44  ; 

45  MISSN1-.; 

46  MISS25-.J 

47  MISS27-.J 

48  MISS28-.; 

49  IF  CASE  <  100  THEN  Z  -  0; 

50  IF  CASE  >-  100  THEN  Z  -  1; 


80 


2 


SAS(R)  LOG  CMS  SAS  5.16 


VM/ CMS  CMS  USER  QSECLSB 


51  IF  Z-l  THEN  CASE-CASE- 100; 

52  DO  I  -  1  TO  33; 

53  X  -  1-4;  XSQR  -  X*X;  XZ  -  X*z;  XSQRZ  -  X*X*Z;  Y  -  WEEK  {i} ; OUTPUT; 

54  END; 


NOTE:  INFILE  X  IS  FILE  CORTICO  DAT  Al 
NOTE:  152  LINES  WERE  READ  FROM  INFILE  X. 

NOTE:  DATA  SET  WORK.TESTC  HAS  5016  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.80  SECONDS  AND  204K. 

55  PROC  CONTENTS; 

NOTE:  THE  PROCEDURE  CONTENTS  USED  0.20  SECONDS  AND  460K  AND  PRINTED  PAGES  1  TO  2. 

56  PROC  PRINTTO  NEW  UNIT-20; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 

57  PROC  SORT  OUT-SCTR; 

58  BY  Z  X  Y; 

NOTE:  DATA  SET  WORK.SCTR  HAS  5016  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  PROCEDURE  SORT  USED  1.13  SECONDS  AND  6924K. 

59  PROC  SUMMARY; 

60  BY  Z  X; 

61  VAR  Y; 

62  OUTPUT  OUT-OVLMN  MEAN-MEAN; 

NOTE:  THE  DATA  SET  WORK.OVLMN  HAS  66  OBSERVATIONS  AND  5  VARIABLES. 

NOTE:  THE  PROCEDURE  SUMMARY  USED  0.82  SECONDS  AND  460K. 

63  DATA  SCORTICO; 

64  SET  SCTR  OVLMN; 

65  BY  Z; 

NOTE:  DATA  SET  WORK. SCORTICO  HAS  5082  OBSERVATIONS  AND  10  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.84  SECONDS  AND  332K. 

66  PROC  PLOT  NOLEGEND  DATA-SCORTICO; 

67  BY  Z; 

68  PLOT  MEAN*X« ' X '  Y*X-'.'  /  VAXIS-0  TO  40  BY  5  OVERLAY; 

69  TITLE  ’CORTICOSTERONE  SCATTER  DIAGRAM'; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.62  SECONDS  AND  460K  AND  PRINTED  PAGES  3  TO  4. 

70  PROC  PRINTTO  NEW  UNIT-21; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 

71  PROC  PLOT  NOLEGEND  DATA-SCORTICO; 

72  PLOT  MEAN*X- ' X '  /  VAXIS-0  TO  40  BY  5; 

73  TITLE  'Mean  Plasma  Corticosterone  Concentrations  Versus  Time'; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.22  SECONDS  AND  460K  AND  PRINTED  PAGE  5. 

74  PROC  PRINTTO  NEW  UNIT-22; 

75  TITLE  'CORTICOSTERONE  ANALYSIS' ; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 
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76  PROC  DATASETS; 

77 

LIST  OF  MEMBERS  BEFORE  UPDATE  OF  DIRECTORY. 


NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

OVLMN 

/DATA 

66 

1 

SCORTICO/DATA 

5082 

1 

SCTR 

/DATA 

5016 

1 

TESTC 

/DATA 

5016 

1 

77 

DELETE  SCTR; 

78 

DELETE  OVLMN; 

79 

DELETE  TESTC; 

LIST  OF 

MEMBERS  AFTER 

UPDATE  OF 

DIRECTORY. 

NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

SCORTICO/DATA 

5082 

1 

NOTE:  THE  PROCEDURE  DATASETS  USED  0.13  SECONDS  AND  460R. 

80  PROC  STEPWISE; 

81  MODEL  Y  -  X  XSQR  Z  XZ  XSQRZ 

82  /  STEPWISE  MAXR; 

NOTE:  THE  PROCEDURE  STEPWISE  USED  0.90  SECONDS  AND  460K  AND  PRINTED  PAGES  6  TO  11. 

83  PROC  PRINTTO  NEW  UNIT-23; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 

84  PROC  REG; 

85  MODEL  Y  -  X  XZ  XSQRZ  /  PARTIAL; 

86  ID  CASE; 

NOTE:  ACOV  AND  SPEC  OPTION  ONLY  VALID  WITH  RAWDATA 

NOTE:  THE  PROCEDURE  REG  USED  2.83  SECONDS  AND  652K  AND  PRINTED  PAGES  12  TO  16. 

87  PROC  PRINTTO  NEW  UNIT-24; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 

88  PROC  GLM; 

89  CLASS  X  Z; 

90  MODEL  Y  -  X  Z  X*Z; 

NOTE:  THE  PROCEDURE  GLM  USED  5. '’5  SECONDS  AND  1036K  AND  PRINTED  PAGES  17  TO  18. 

91  PROC  PRINTTO  NEW  UNIT-25; 

92  * - * 

93  *  * 

94  *  to  o  in  the  tables  listing  the  variance  inflation  * 

95  *  fact,..*,  influence  statistics,  and  tolerances,  the  * 

96  *  following  SAS  statements  were  used  in  this  partition:  * 

97  *  * 

98  *  PROC  REG;  * 


99  *  MODEL  Y  -  X  XZ  XSQRZ  /  TOL  VIF  INFLUENCE;  * 

100  *  ID  CASE;  * 

101  *  OUTPUT  OUT-RCORTICO  P-PREDICT  R-RESID  STUDENT-STUDENT  * 

102  *  * 

103  * - *; 


NOTE:  THE  PROCEDURE  PRINTTO  USED  0.04  SECONDS  AND  332K. 

104  PROC  REG; 

105  MODEL  Y  -  X  XZ  XSQRZ  /  I  SSI  SS2  STB  COVB  CORRB  SEQB  COLLIN 
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106  COLLINOINT  ACOV  P  R  CLM; 

107  ID  CASE; 

108  OUTPUT  OUT=RCORTICO  P-PREDICT  R-RESID  STUDENT-STUDENT; 

NOTE:  THE  DATA  SET  WORK.RCORTICO  HAS  5082  OBSERVATIONS  AND  13  VARIABLES. 

NOTE:  THE  PROCEDURE  REG  USED  11.63  SECONDS  AND  652K  AND  PRINTED  PAGES  19  TO  126. 

109  PROC  PRINTTO  NEW  UNIT-26 ; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 

110  PROC  PLOT  DATA-RCORTICO; 

111  PLOT  RES ID*X-'*'  ; 

112  PLOT  RESID*PREDICT- ' * ' ; 

113  PLOT  STUDENT*X- ' * ' ; 

114  PLOT  STUDENT*PREDICT- ' * ' ; 

115  TITLE  'CORTICOSTERONE  ANALYSIS' ; 

NOTE:  THE  PROCEDURE  PLOT  USED  2.21  SECONDS  AND  460K  AND  PRINTED  PAGES  127  TO  130. 

116  PROC  PRINTTO  NEW  UNIT-27; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  332K. 

117  PROC  PLOT  DATA-RCORTICO; 

118  BY  z; 

119  PLOT  RESID*CASE- 1  * 1  /  HAXIS-1  TO  90  BY  5; 

120  PLOT  STUDENT*CA3E« ’ * ’  /  HAXIS-1  TO  90  BY  5; 

121  TITLE  ’CORTICOSTERONE  ANALYSIS' ; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.79  SECONDS  AND  460K  AND  PRINTED  PAGES  131  TO  134. 
NOTE:  SAS  USED  6924K  MEMORY. 

NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 

PO  BOX  8000 

CARY,  N.C.  27511-8000 


APPEMDIX  I 


STEPWISE  AND  MAXIMUM  R2  REGRESSION  PROCEDURES 
USED  TO  BUILD  CORTICOSTERONE  MODEL 
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GENERAL  LINEAR  MODELS  PROCEDURE 
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this  term  contains  both  sum-of-squares  pure  error 
and  sum-of-squares  lack-of-fit. 
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